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We introduce a generic approach to study interaction effects in diffusive or chaotic quantum dots 
in the Coulomb blockade regime. The randomness of the single-particle wave functions induces 
randomness in the two-body interaction matrix elements. We classify the possible induced two- 
body ensembles, both in the presence and absence of spin degrees of freedom. The ensembles 
depend on the underlying space-time symmetries as well as on features of the two-body interaction. 
Confining ourselves to spinless electrons, we then use the Hartree-Fock (HF) approximation to 
calculate HF single-particle energies and HF wave functions for many realizations of the ensemble. 
We study the statistical properties of the resulting one-body HF ensemble for a fixed number of 
electrons. In particular, we determine the statistics of the interaction matrix elements in the HF 
basis, of the HF single-particle energies (including the HF gap between the last occupied and the 
first empty HF level), and of the HF single-particle wave functions. We also study the addition of 
electrons, and in particular the distribution of the distance between successive conductance peaks 
and of the conductance peak heights. 



PACS numbers: 73.23.Hk, 05.45.Mt, 73.63.Kv, 73.23.-b 



I. INTRODUCTION 

In this paper, we explore a generic approach to inter- 
action effects in diffusive and/or chaotic quantum dots 
in the Coulomb blockade regime. To set the stage, we 
describe in this Introduction some previous work, and 
discuss the motivation of our own investigation. 

Electron transport in quantum dots that are strongly 
coupled to leads (open dots) has been successfully de- 
scribed within a single-particle framework. Electron- 
electron interactions affect the dephasing rate of the elec- 
trons at finite temperature, but otherwise electrons in the 
vicinity of the Fermi energy can be described as non- 
interacting quasi-particles. The single-particle dynam- 
ics in the dot can be affected by disorder or by the dot's 
boundaries. In a disordered dot, the electron moves diffu- 
sively. If the disorder is weak, the statistical fluctuations 
of the single-electron spectrum and wave functions can 
be described by random-matrix theory (RMT).-^- Sim- 
ilarly, in a ballistic dot with irregular boundaries, the 
classical dynamics are mostly chaotic and the statistical 
quantal fluctuations also follow RMT. In such diffusive 
and/or chaotic open dots, the conductance exhibits "ran- 
dom" but reproducible fluctuations versus, e.g., gate volt- 
age or magnetic field that can be successfully described 
by RMTiS 

The situation changes as the coupling of the dot to the 
leads is reduced and the dot becomes almost isolated. 
The charge on the dot becomes quantized, the conduc- 
tance as a function of gate voltage displays sharp peaks 



(the Coulomb blockade resonances), and the Coulomb in- 
teraction between electrons becomes important.^ In this 
Coulomb blockade regime, the simplest model for a quan- 
tum dot is the constant interaction (CI) model. The 
interaction is modeled by the classical charging energy 
of a system with capacitance C. For a fixed number of 
electrons n on the dot, the interaction of this model is 
constant, and the model essentially reduces to a single- 
particle model. While certain statistical properties of 
the conductance peak heights can be described within 
the CI plus RMT model)^^ other observables, most no- 
tably the distribution of the distance between successive 
conductance peaks (a distance known as the peak spac- 
ing), deviate significantly from the predictions of the CI 
plus RMT model pMilili These deviations have been un- 
derstood, at least qualitatively, to follow from residual 
interaction effects beyond the charging energy-'^J-P Such 
residual interaction effects are at the center of the present 
investigation. 

In a weakly diffusive and/or chaotic dot, the "random- 
ness" of the single-particle eigenfunctions induces ran- 
domness of the interaction matrix elements. The sta- 
tistical fluctuations of these interaction matrix elements 
(evaluated in the non-interacting basis) were calculated 
for a dot with a large Thouless conductance g^^ and 
were found to be suppressed. In the limit g -^ oo, only 
a few interaction terms survive, leading to the so-called 
universal Hamiltonian.^^'^- For spinless electrons, the in- 
teraction part of the universal Hamiltonian is composed 
of just the charging energy term. This shows that the CI 



model is appropriate for spinless electrons and in the limit 
g —f oo. When the electron spin is taken into account, 
there is also an exchange interaction term, and, in the 
absence of a time-reversal-symmetry breaking magnetic 
field, a Cooper-channel-like term. The exchange inter- 
action term was found to have important effects on the 
finite-temperature peak spacing distribution. ^^-^^ For fi- 
nite g, further interaction terms exist and can be calcu- 
lated within a systematic expansion of the residual inter- 
action in inverse powers of g. Using a two-body effective 
screened interaction in the limit of a small gas constant 
rs , one can show that for a diffusive dot the average and 
the standard deviation of these residual interaction terms 
are at least of the order A/g, where A is the single- 
particle mean level spacing, while for a chaotic dot the 
standard deviation of the residual interaction terms is 
of the order A^/lng/g/i^t^ For finite values of g, such 
terms affect the ground-state energy of the dot and var- 
ious observables. In particular, the position of the con- 
ductance peak height at low temperatures is determined 
by the change of the ground-state energy upon the ad- 
dition of an electron to the dot (addition energy). Thus, 
finite-(7-corrections to the universal Hamiltonian affect 
the peak-spacing distributionr-^^J^ 



Most studies of interaction effects have been limited 
to specific models of a dot. These models include an 
Anderson model for a diffusive dot or a billiard model 
for a ballistic chaotic dot, plus Coulomb or Hubbard-like 
(short-range) interactions. For a small number of elec- 
trons, exact numerical solutions are possible.^'^'' For a 
larger number of electrons, such models have been stud- 
ied in the Hartree-Fock (HF) approximationi2ii2iS£ and 
in a density functional approach. ^^'^^-^'^ These models 
contributed significantly to our understanding of interac- 
tion effects in almost isolated dots, but they also contain 
non-generic features. 

Our approach to interaction effects in diffusive and/or 
chaotic, almost isolated quantum dots aims at a generic 
understanding of interaction effects beyond the g —^ oo 
limit. Our starting point is a Hamiltonian consisting of 
a one-body and a two-body part. The one-body part 
describes chaos and/or disorder and is, within an energy 
interval of width « ^A, described by the appropriate 
random-matrix ensemble, i.e., the orthogonal ensemble 
for conserved time-reversal symmetry and the unitary 
ensemble for broken time-reversal symmetry. The two- 
body part is given by the Coulomb interaction. However, 
the randomness of the single-particle wave functions in- 
duces a randomness of the two-body interaction matrix 
elements when the latter are written in terms of the eigen- 
functions of the random one-body part of the Hamilto- 
nian. This results in a generic ensemble of the two-body 
interaction that we derive and use in our studies. 

The paper consists of two parts. In the first part, we 
derive the statistical properties of the two-body matrix 
elements (the "induced two-body ensemble"). We find 
that these properties are determined both by the under- 
lying space-time symmetries of the system and by fea- 



tures of the two-body interaction. We classify the pos- 
sible forms of the first two cumulants of the interaction 
matrix elements. This is accomplished by writing these 
cumulants in a covariant form, ensuring invariance under 
a change of the single-particle basis. The first moments 
(average values) reproduce the universal Hamiltonian in 
the limit 5 — > oo. The second cumulants are charac- 
terized by a constant m^ {u^ cx ^ I g^ for a diffusive 
dot and m^ oc A^ In gj g^ for a chaotic dot) . The invari- 
ance classes of the second cumulants are determined by 
the symmetry properties of the interaction matrix ele- 
ments under permutations of the single-particle orbitals. 
These symmetries are different for, e.g., a contact in- 
teraction or a short- but finite-range interaction. We 
note that the two-body part of one class of the ensem- 
bles which we construct, coincides with the two-body 
embedded Gaussian ensemblesi^ These ensembles were 
introduced in the nuclear shell model to study the va- 
lidity of RMT for the statistical description of nuclear 
spectrai2S*2& The spectral properties of these embedded 
ensembles were recently analyzed using an eigenvector 
expansion of the second moments^iS^ Combined with 
a random-matrix one-body part, such a two-body en- 
semble was used to study interaction effectsSSiffi and 
ground-state magnetization^ in quantum dots with a 
small number of electrons. The spin structure of a system 
with a spin-conserving random interaction was studied 
in Ref . [S^. 

In the presence of spin degrees of freedom, one needs 
to use non-antisymmetrized two-body matrix elements, 
and our classification of the second cumulants is gener- 
ally done for such elements. The spinless case requires 
only antisymmetrized matrix elements. Their second cu- 
mulants are obtained directly from the cumulants of the 
non-antisymmetrized elements. The two-body ensem- 
bles that correspond to a contact interaction exist only 
in the presence of spin. 

In the second part of the paper, we use the Hartree- 
Fock (HF) approximation to work out interaction effects 
for some of the ensembles introduced in the first part. It 
is difficult to study the properties of ensembles involving 
a two-body interaction. Exact numerical solutions can 
be obtained only for dots with a small number of elec- 
trons. For a larger number of electrons, it is necessary 
to make approximations. The HF approximation serves 
this purpose. It has the advantage of being based on a 
single-particle picture (optimized to take into account in- 
teraction effects), and one can still use concepts familiar 
from the independent-particle approach. While the HF 
approximation has been studied in the context of specific 
models for an interacting dotji2ii2i2S the present work is 
free of any model-dependent effects and constitutes a 
generic statistical HF approach. 

Our HF studies are limited to spinless electrons. We 
investigate, in particular, the statistical properties of the 
interaction matrix elements in the basis of HF single- 
particle eigenfunctions, and of the HF single-particle 
eigenvalues and eigenfunctions. The interest in such work 



derives, among others, from generic studies of quantum 
dots with a large number of electrons. These often use 
a HF single-particle basis (rather than a basis defined 
in terms of non-interacting electrons). ^^ It is then of- 
ten assumed that the statistical properties of the matrix 
elements in the HF basis are similar to those of the non- 
interacting basis. ^ Our work permits us to test such as- 
sumptions, whereas previous studies (done in the frame- 
work of specific models of a dot with interactions) have 
focused on the spectral statistics of the HF levels only. 
We find that, in the limit of small interaction strength u, 
the average interaction matrix elements in the HF basis 
acquire a v? correction (this explains the increase of the 
average HF gap with u), while the second cumulants re- 
main very close to their input values. As for the spectral 
statistics, our generic studies confirm the conclusions of 
previous studies. The nearest-neighbor spacings of filled 
and empty HF levels follow the Wigner distribution of 
RMT, while the distribution of the HF gap deviates from 
this distribution. An open problem has been whether 
correlations of filled and empty levels follow or deviate 
from RMT. Our studies of an appropriate spacing corre- 
lator indicate that the gap weakens the correlations be- 
tween spacings on its opposite sides. On the other hand, 
this spacing correlator is found to follow RMT within 
the filled or empty levels alone, adding further support 
to the conjecture that the filled and empty levels sepa- 
rately satisfy RMT statistics. We also show that the HF 
ensemble satisfies orthogonal (unitary) invariance, and 
therefore its wave function statistics is that of RMT. 

We use the statistical HF approach to study the statis- 
tics of observables, and in particular of the peak-spacing 
and peak-height distributions. The results can be easily 
interpreted in Koopmans' limit j2^ which assumes that the 
HF wave functions remain unchanged upon the addition 
(or removal) of an electron. Koopmans' limit has been in- 
vestigated in previous HF studies jiSiiSiSS but the present 
approach again has the advantage of being generic. It al- 
lows us to define the limitations of Koopmans' approach, 
and to go beyond that limit. In particular, we find that 
the peak-spacing distribution can be described as a con- 
volution of a Wigner distribution (with an appropriate 
mean level spacing) and a Gaussian whose width is ap- 
proximately the standard deviation of a diagonal interac- 
tion matrix element. On the other hand, we find that the 
peak-height distribution is not affected by the residual 
interaction. 

The outline of this paper is as follows: In Section^ we 
discuss the induced two-body ensembles. We classify the 
possible ensembles for the orthogonal and unitary sym- 
metries in Sections III B II and III B 21 respectively, both in 
the presence and absence of the spin degrees of freedom. 
The possible ensembles are related to different types of 
two-body interactions, as characterized by the symme- 
tries of their matrix elements under permutations of the 
orbitals. In Section Fill Al we discuss the HF approach to 
disordered (or chaotic) systems in the framework of the 
many-body ensembles. A perturbative approach that 



provides a closed expression for the average HF matrix 
elements for a fixed single-particle spectrum is described 
in Section Fill Bl Numerical studies of both the first and 
second moments of the interaction matrix elements in the 
HF basis are presented in Section IIII CI The statistics 
of the HF single-particle energies including filled levels, 
empty levels and the gap are studied in Section IIIIDI 
while the statistics of the HF wave function components 
are studied in Section IIII El Finally in Section IIVI we 
discuss the addition of electrons in the HF approach to 
disordered systems, and in particular the peak-spacing 
and peak-height distributions. 



II. DISORDERED SYSTEMS WITH 
INTERACTIONS 

We consider a disordered system of interacting elec- 
trons, e.g., a diffusive quantum dot. For simplicity, we 
first discuss the case of spinless electrons. The spin de- 
grees of freedom are introduced later in this Section. For 
a given realization of disorder, we describe the system 
by a Hamiltonian with a one-body part and a two-body 
interaction. 
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The m orthonormal single-particle states \i) — a]|— ) 
(i — 1, . . . , m) with |— ) the vacuum state, form a fixed 
basis with m '3> 1- The one-body part h^-' of H 
stands for kinetic energy plus disorder potential, while 
Vij-ki represents the (screened) Coulomb interaction. The 
Vij-ki = (ijlvlkl) are non-antisymmetrized matrix ele- 
ments of the two-body interaction Qij) is a product state 
of particle 1 in state \i) and particle 2 in state \j)). For 
spinless fermions, we can rewrite the interaction in terms 
of antisymmetrized matrix elements w^.^.; = Vij-ki — Vij-ik 
as 
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Disorder is taken into account by postulating that h 
is random. The way this randomness is modeled turns 
out to be rather immaterial. Indeed, we recall that for 
weak disorder and in the long- wavelength limit, all mod- 
els for non-interacting disordered systems yield the same 
supersymmetric field theory, the non-linear supersym- 
metric sigma modeL34 This field theory contains one es- 
sential parameter, the dimcnsionless Thouless conduc- 
tance g. We will work out properties of the system de- 
scribed by the Hamiltonian of Eq. |5J| in terms of an 
expansion in inverse powers of g. We will carry the ex- 
pansion up to and including terms of order l/g^- 

The non-linear sigma model shows that to zeroth order 
in 1/(7, h-, is equivalent to one of the canonical ensembles 
of RMT, the Gaussian orthogonal (GOE), unitary (CUE) 



or symplectic (GSE) ensemble of random matrices, de- 
pending upon the symmetry of the problemi^ Correc- 
tions calculated by using the expansion in inverse powers 
of g yield deviations from predictions of RMT but do not 
violate the fundamental symmetries of the problem (in- 
variance under orthogonal, unitary or symplectic trans- 
formations of the basis of single-particle states). Indeed, 
the derivation of the non-linear sigma model is based 
upon these symmetries, which are thus deeply embed- 
ded in the model and are valid beyond canonical RMT. 
This fact will be important as we will heavily rely on 
these symmetries to determine the effective form of the 
Hamiltonian H. We will focus on the cases of orthogonal 
and unitary symmetry and will not consider symplectic 
ensembles in the sequel. 

The randomness of /i^- induces a degree of random- 
ness into the two-body matrix elements Vij-ki- To see 
this we recall that every observable O is defined as an 
average over the ensemble /i^ . Because of the invari- 
ance of that ensemble under orthogonal (or unitary or 
symplectic) transformations of the single-particle basis, 
O must likewise be invariant under all such transforma- 
tions. Therefore, O can depend on the two-body interac- 
tion Vij-ki only through invariants constructed from the 
matrix elements Vij-ki- Other properties of the Coulomb 
interaction that are not encapsulated in the invariants 
are not relevant. On a formal level, this fact allows us 
to define equivalence classes of two-body interactions. 
Members of the same class share the same set of invari- 
ants although the interactions may differ in other proper- 
ties. All two-body interactions which belong to the same 
class as the Coulomb interaction will give rise to identical 
average properties of the ensemble ^ . These members 
may be said to form an induced ensemble of two-body 
interactions. The word "induced" indicates that the en- 
semble owes its existence to the randomness of the one- 
body part /i- ■ . Thus, we deal effectively with an induced 
two-body random ensemble although the Coulomb inter- 
action per se is fixed and not random at all. We aim at 
reformulating the ensemble of Hamiltonians H in Eq. (Q 
in such a way that only the relevant features (i.e., the 
invariants) of Vij-ki are kept. In doing so, we do not fol- 
low strictly the line of argument just given because we 
perform separately the averages over the eigenf unctions 
and eigenvalues of /i^ ■ . The induced two-body ensemble 
will result from the average over the eigenfunctions. 

We accordingly transform H in Eq. Q to the basis \a) 
of single-particle eigenstates of ft.^"-' . This yields 



H = y^£aa]^Qa 
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In an interval of length ~ gA (where A is the aver- 
age single-particle level spacing), the eigenvalues e^ obey 
canonical random-matrix statistics. Because of the ran- 
domness of the single-particle eigenfunctions, the ele- 
ments v^a.^g are also random. Moreover, the e^'s are 
uncorrelated with the eigenfunctions and, thus, with the 



Ua/3;75's. This last property holds strictly if h^*^^ stands 
for the GOE or CUE, i.e., in the limit g ^ oo. For a 
diffusive system, it holds up to and including terms of 
order 1/g'^. 

The number of single-particle levels a is usually cho- 
sen to be finite. Thus, the interaction in Eq. (Q is gen- 
erally an effective interaction, e.g., a screened Coulomb 
interaction. Such an effective interaction was, for exam- 
ple, derived in Refs. lollla within an interval of ~ g levels 
around the Fermi energy of the dot. The authors used the 
random-phase approximation (RPA) in the limit where 
the gas constant Ts <C 1. The effective interaction is 
then given by the unscreened zero-momentum Fourier 
component of the Coulomb interaction plus a short-range 
screened Coulomb interaction (whose exact form depends 
on the dimensionality and geometry of the system). In 
general, for a finite system this effective interaction also 
includes surface charge termsiS These can be eliminated 
by assuming, for instance, periodic boundary conditions 
and will not be taken into account in the sequel. 

If the spin of the electrons is taken into account, the 
Hamiltonian ^ is replaced by 
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where a = ± describes a spin up/down electron and \i) 
denotes an orbital state. The one-body Hamiltonian h'^^^ 
is spin-independent (i.e., there is no spin-orbit coupling). 
The interaction matrix elements are also assumed to be 
spin-independent as is the case for the Coulomb interac- 
tion. 

In analogy to the transformation leading from Eq. (Q 
to Eq. 0, the ensemble Q can be written in terms of 
the orbital eigenfunctions \a) of the one-body Hamilto- 
nian /i^"'. Since the transformation from the basis of 
states \i) to that of the states |a) does not involve the 
spin degree of freedom, the single-particle energies e^ are 
spin-degenerate, and the interaction matrix elements re- 
main spin-independent. Thus, when the spin is included, 
Eq. (j3Jl is replaced by 



i/ = ^ taal^aa^a + ^ E ^"/^i 
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In contrast to the spinless case where H could be ex- 
pressed in terms of the antisymmetrized matrix elements 
V , Eq. (O contains the full non-antisymmetrized inter- 
action matrix elements. 

We now work out the first and second moments of 
Vai3--y5 by averaging over the eigenfunctions |a). 



A. Mean Values 

We first calculate the mean value Va/s-^js of the two- 
body matrix elements, invoking orthogonal or unitary 



invariance. In the orthogonal case there are three in- 
variants: dofySps, Sas^fSf, and SafsSjs- In the unitary 
case, only the first two invariants exist. We consider 
first the case of antisymmetrized matrix elements and 
include spin below. For the antisymmetrized element 
'^aB-jS = Vaf3--yS — Vafj-S'y, the first two invariants combine 
to SajSps — SasSffj, while the third one gives no contribu- 
tion. We conclude (for both the orthogonal and unitary 
symmetries) that 

""afS-jS = VoiSa-fS/SS ~ SasS/S^y) ■ (6) 

The average interaction in |j3l is then given by 



'^ = 2""°^^ ""'-' ' 
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where n is the particle-number operator. For a fixed 
number n of fermions, the average interaction is sim- 
ply a constant voin^ — 'n)/2. For a quantum dot, the 
zero-momentum Fourier component of the Coulomb in- 
teraction is unscreened and its contribution to vq is the 
charging energy e^/C where C is the dot's capacitance. 
The contribution of the screened short-range Coulomb 
interaction can be estimated using the diagrammatic ex- 
pansion and is of order A. For finite g, there are correc- 
tions of order \/g. 

For the case with spin, all three (two) orthogonal (uni- 
tary) invariants now contribute to the value of the aver- 
age interaction. 



ya/3;7<5 — V^SaySpS + JsSasSpj + JcSasSf3j 



(8) 



Here wq, Js and Jc are constants and the last term on 
the r.h.s. is absent for the unitary case. Using Eq. (jSJ in 
Eq. |SJl we obtain an average interaction of the form 

v^\{vo- Js/2)n2 - (wo/2 - Js)n - J,s' + J.f^f . (9) 

In Eq. linji, S = i J2a J2cra' O-ta^^cra' O-aa' is the total Spin 

operator of the dot, where a is the vector of the three 2x2 
Pauh matrices. The operator T^ = J2a ^I+'^I- creates 
coherent pairs of spin up-down electrons. Adding the av- 
erage interaction J^ to the one-body part of Eq. jSJ and 
taking the limit g — > oo, we obtain the universal Hamilto- 
nian of the quantum dot, with ■^o = Vap-ap, Js = Vap^pa, 
and Jc = Vaa;pp the direct, exchange and the Cooper 
channel interaction strengths, respectivelyi ^^i^'^ We note, 
however, that Eq. ^ is valid beyond this limit. For finite 
values of g, corrections to the various constants can be 
obtained by the diagrammatic expansion in 1/g. Unlike 
the charging energy term (which for a fixed number of 
electrons is simply a constant), the exchange and Cooper 
channel interactions are represented by non-trivial oper- 
ators. 



B. Second Moments 

The second moments of the interaction matrix ele- 
ments depend upon the symmetry of the ensemble, and 



upon the generic form of the two-body interaction. It is 
easier to calculate first the second moments of the non- 
antisymmetrized matrix elements Vap-^s for the two- 
body ensemble in the presence of spin degrees of free- 
dom. The spinlcss case is then obtained by considering 
the corresponding moments of the antisymmetrized ma- 
trix elements v^p.^g. 

To simplify the notation, we assume in the following 
that the average interaction has been separated out and 

put Vap-jS = 0. 



1. Orthogonal induced two~body ensembles 

In the case of orthogonal symmetry the two-body in- 
teraction matrix v^p-rys is real and symmetric, Vap-^s = 
'"ap-fS ~ '^■yS;ap- Depending on the form of the interac- 
tion, we distinguish three induced ensembles. 
(i) For a schematic effective interaction that is a S- 
function (i.e., contact interaction), the matrix element 
Vap-yS = VAJ drilja{'r)ipp{r)tpt{'r)ips{'r) {V is the vol- 
ume and A is the single-particle mean-level spacing) is 
invariant under all 24 permutations of the four indices 
aP"fS. An ensemble of two-body matrix elements that 
satisfies the orthogonal invariance condition and is con- 
sistent with the symmetries of Vap-^-yS, is expected to have 
to leading order in 1/g (see below) a second moment of 
the formJ^ 



Vap;ySV^j.v;pa — 3u {5ai_i5 p^SypSsa 
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where . . . stands for the sum of terms obtained by 
all 23 permutations of (xvpcr. Non-linear sigma model 
calculations fii as well as diagrammatic calculations in 
weakly disordered systems using a screened Coulomb in- 
teraction, show that u^ (x A^/g^ (see below). We have 
written down only invariants that arise by pairing in- 
dices appearing in different matrix elements. Invariance 
requirements alone would also allow terms involving Kro- 
necker symbols like Sap etc. However, in diffusive systems 
such invariant terms are at least of third order in 1/g and 
are, therefore, ignored. An analogous statement applies 
to all the other induced two-body ensembles listed be- 
low, both for the orthogonal and the unitary cases, and 
will not be repeated there. 

For a diffusive dot, the quantity u^ in Eq. (|10|l can be 
calculated in the diagrammatic approach. If the energy 
difference between any pair of levels from the set a(3^5 is 
smaller than the Thouless energy gA, Eq. (|10|l is found 
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fdrfdv'n\r,r')/V^ = c^, (11 
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independently of the energy differences of 
the single-particle levels. Here H(ri,r2) — 

(l/A/^)^Q_,o0Q(ri)(/.Q(r2)/i?Q2 is the diffusion 
propagator in the finite dot, with 0q the eigenfunction 
of the diffusion operator corresponding to the eigenvalue 



_DQ^. The proportionality constant c in Eq. I|ll() 
depends on the geometry of the dot. For a cubic dot of 
length i in d dimensions, c = ;^ X^n^^o l/l'^l^ where n is 
a vector of d non-negative integers (and the dimension- 
less conductance g is defined by gA — 2'KhD / L?') .^^ For 
a 2D circular dot of radius E, c = 2[2, ,„ x^XV''^ « 0.67 
where xi_m are the zeros of the derivative of the Bessel 
function of order I (and the dimensionless conductance is 
defined by gA — 27ThD/R'^)i^ In the general case where 
some of the energy differences between pairs of levels 
are greater than gA, the coefficient c in Eq. (|11|) is no 
longer constant but depends on these energy differences. 
Eq. l(TT|l can be extended to a balHstic dot. This is 
done with the help of a supersymmetric sigma model, 
obtained through the addition of weak disorder with fi- 
nite correlation length. ^^^ In the ballistic case^ 



c'^Hc"g) 



where g = 7r(n/2)^/^ is the ballistic Thouless conduc- 
tance, n is the number of electrons, and c' = 3/4 is 
a geometry-independent constant. The constant c" de- 
pends on the geometry of the dot. For a circular ballistic 
dot, c" « 0.81. 

From Eq. (|10|l . the variances of off-diagonal and di- 
agonal matrix elements are given by {a,P,"f,6 are all 
different) 
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For spinless fermions, we have v^^ g = Vai3--yS — 
Vap-Sj = because of the symmetries satisfied by the u's. 
Therefore, there is no non-trivial ensemble of antisym- 
metrized matrix elements that corresponds to Eq. H10|l . 
and the ensemble IjlOl) exists only when the spin degrees 
of freedom are taken into account. 



(ii) For a generic two-body local symmetric interaction w(ri, r2) = w(r2j ^2)^ the matrix elements satisfy the relations 



The second moment of the ensemble is expected to have the form 

^ap^l3nu^1fj.^Scr + ^av^ PpSfa^Smu + ^ap^PaS-fpSsv + ^aa^ fifi^^u^Sp) ; 



(14) 



(15) 



where the eight terms correspond to the eight permutations of fivpa that leave v^^p„ invariant (see Eqs. H14|) '). The 
form (|15|l is also obtained in the diagrammatic approach to order 1/ g^ when an RPA two-body screened interaction 
w(ri,r2) is used. In the limit r^ ^ 1, some of the diagrams that do contribute for a (5-function interaction, give a 
negligible contribution for a local finite-range interactiouj^ this leading to Eq. H15|l instead of Eq. H10|l . The factor 
v? is still given approximately by Eq. (|ll|l . but now 

<y'^[Vap--is) =U^ , 0-^(Wa/3;Q/3) = 4u^ , Cr^ (Wq^;^q) = 2u^ , CT^ (WaQ;aa) = 8u^ . (16) 

The direct {vap-^-ys) and exchange {vai3;S-y) matrix elements are uncorrelated. The ensemble of antisymmetrized matrix 
elements v^p.^^ follows from Eq. (fT5|l using w^^..,xW 
We find 



af3;f5 pi^'tPry 



j^ 1^ 



In particular 



'"ap-'y5'^py\pcr — 2u {SapSp^S-^pSsa ~ SapSp^S-fcrSsp + Sa„Sf3^S^^Ssp — 5a^Sf3^S^pSsa- 

+SapSf3aSjpSsu ~ SaaSppS^f^Ss,, + SaaSppSj^Ssp — SapSf3c,5ji,Ssp) 

+U {SapS/3^6jfj,6scr — SaaSpi,S^pSsp + Sau^Pp^ja^Sp — SaijSfjcrSjpSsp 

+SapSpaSjp6si, — SapSppS^a^Su + SaaSppS^^Ssp — 6ap6pp,5j^5sa 

+5aa5pv5'yp5&p — SapSp^SjcrSsp + 5ay5pa^^p^&p — SaijSppSj^Ssa 

+SapSppSj^5scr — SapS/3crS^i,Ssp + SapSppS^crSsu ~ ^aa^ Pp^-ip^&v) ■ 



^''k^tp-a^) = 2^2 , a^{v%,^i,) = ^v? 



(17) 



(18) 



(iii) The matrix elements of a generic non-local two-body interaction u(ri, r2; r'^^, ra) = w(r2, ri; r2, r'^^) (we define 
w(ri,r2;ri,r2) = (ri,r2|-y|ri,r2)) satisfy 



VaP-^-ii = ViSa-S-y = VjS-ap = VSj-^po 



(19) 



Therefore, we expect for the ensemble average 

i.e., the two-body Gaussian orthogonal ensemble. The corresponding variances are given by 

Cr^(Wa/3;75) = W^ , Cr^ (wQ/3;a/3) = 2u^ , CT^ (Wc/3;/3a) = 2u^ , <J^ {Vaa;aa) = 4:u'^ ■ (21) 

For the antisymmetrized interaction matrix elements we expect accordingly 



^a/3;7<5^M!^;pT — ^U {SafiSp^S^pSsa — 5ap5 piyS-yahp + S^u^ Pt^S^aSsp ~ S^vS ppS^ pSsa 

+ 5ap5po5-ip5su ~ 5aa5j3p5^p,5sv + ^aa^ PpS^v^S^ — S^pS po^^u^Sp) ■ (22) 

In particular 

a'{v^f,.^^s)^2u\ a\vtp,^p) ^ An" . (23) 

These relations should be compared with Eqs. p8(l . It is not clear whether there is an interaction model for which 
disorder averaging leads to the ensemble H22(l . although this ensemble seems to be generic for non-local interactions. 

2. Unitary induced two-body ensembles 

For the unitary symmetry the two-body interaction matrix is complex Hermitcan, Vap.^^s — v*s-ai3- There are two 
possible two-body ensembles: 
(i) The matrix elements Vap.^g satisfy the relations 

The ensemble consistent with relations (|24|) is 



V*ap:'^S'"p,i^\p(y = 2m {SapSpuS-^pSsa + Sai^SppS^pSgc + 6apSpi,SjaSsp + 5au5pp5^„5sp) ■ (25) 

In the disorder basis such an ensemble is realized for a J-function interaction with u^ given by Eq. (|ll|l . Variances of 
diagonal and off-diagonal elements are given by 

0'^(Wq/3;75) = 2m^ , a"^ {VaP;ap) ^ 2u'^ , a^ (Vap-pa) ^ 2u'^ , Ct"^ (Vaaiaa) = Su^ . (26) 

Relations (|24ll lead to v^a g — 0, so there is no non-trivial ensemble of antisymmetrized matrix elements that 

corresponds to Eq. H25(l . 

(ii) The matrix elements satisfy relations that are typical for a symmetric two-body interaction (local or non-local) 

VaP;jS = Vpa-S^ = V*^s.^i3 = w|^.^„ . (27) 

A corresponding induced two-body random matrix ensemble obeys 



^a/3;75"M"^;P'7 = "" {^ap5p^5^p5sa + 5au5 ppS^^Ssp) , (28) 

i.e., the two-body Gaussian unitary ensemble. Such an ensemble is realized in the disorder basis by an RPA two- 
body screened interaction in the limit r^ <C 1, with v? given approximately by Eq. pi|l . The variances of the matrix 
elements are 

<y'^{Vap-ns) = U^ , 0-^(Wa/3;a/3) = U^ O-^iVap-pa) = "^ , Cr"^ (Vaa^aa) = 2u^ . (29) 

For the antisymmetrized interaction we find 



Vap-^^S'Vpu-.pa = 2m [SapSpiyS-^pSsa - SapSp^S^a^Sp + 5a,^5 ppS-^cr^Sp - SaiySppS-fpSsa) , (30) 

and in particular 

<j'iv^p,,s)^2u\ a\vtp.^^p)^2u' . (31) 
I 

In the orthogonal case there was a third ensemble, but in the unitary case this third ensemble coincides with the 
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second one because the unitary symmetry relations l|27() 
are common to both local and non-local interactions. 

We reiterate that in the spinless case, the only non- 
trivial induced ensembles are (ii) and (iii) in the orthog- 
onal case (Section IIIB 1|) . and (ii) in the unitary case 
(Section HIEIJ. 



in the numerical studies of the statistics of HF matrix 
elements reported below, such correlated terms do not 
occur, and we use the fully random two-body ensemble. 



C. Induced versus true two— body ensembles 

To simplify the discussion, we disregard the average 
part of the interaction in the ensembles © and ||SJ) and 

put Val3--i& = 0. 

So far, we have calculated only the second moments 
of the two-body interaction. What about higher mo- 
ments? Inspection of the relevant diagrams shows that 
all higher-order cumulants contribute terms of order \/ g^ 
with fc > 3 and are, thus, negligible in the present con- 
text. This does not necessarily imply that the induced 
two-body random ensembles are Gaussian. (It is con- 
ceivable that the kth. cumulants are systematically of or- 
der 1/g^- Rescaling the matrix elements by the factor 
g would then yield a variable which is not Gaussian). 
On the other hand, numerical results in weakly diffusive 
systems do suggest that the interaction matrix elements 
have Gaussian distributions ii^ This can be attributed to 
a central-limit theorem: each matrix element is an inte- 
gral over a product of four wave functions, and each wave 
function is a continuous random variable. 

Our results might suggest that instead of the induced 
random two-body ensembles studied so far, we may use 
true random two-body ensembles. In the latter, the 
matrix elements Vij-^ki are Gaussian-distributed random 
variables with mean value zero and second moments 
given by Eqs. (HUJl, ^^, (EOl, (E3 or (EHIl, as the case 
may be. These relations apply in any single-particle ba- 
sis. In contrast to the induced random two-body ensem- 
bles, the statistical properties of the true two-body ran- 
dom ensembles are inherent rather than induced by the 
randomness of the one-body part of H . Are these two 
sets of ensembles pairwise equivalent? Are we allowed 
to replace a Hamiltonian with a random one-body and a 
fixed two-body interaction by a generic Hamiltonian con- 
taining both a random one-body and the random two- 
body interaction just defined? The answer is no, even 
though both Hamiltonians yield the same distributions 
for the two-body matrix elements. This is because in the 
induced two-body ensembles, randomness is due to the 
eigenfunctions \a) of N-^'. Therefore, the two-body ma- 
trix elements Vap-^-yS are correlated with all other expres- 
sions that depend on the |q;)'s while the matrix elements 
of the true two-body random ensemble arc independent 
random variables that are not correlated with other parts 
of the Hamiltonian. Whenever an observable contains in 
addition to the Ua/3;7(5's other expressions which depend 
upon the |a)'s, the ensemble averages of this observable 
taken over the induced two-body ensemble and over the 
fully random two-body ensemble, will differ. However, 



III. THE HARTREE-FOCK APPROACH TO 

DISORDERED SYSTEMS: FIXED NUMBER OF 

ELECTRONS 



The theoretical treatment of a disordered system with 
interactions poses severe difficulties. This is why we use 
the HF approximation. We express observables in terms 
of the HF single-particle eigenvalues and eigenfunctions. 
We determine the latter by using the HF approximation 
for every realization of disorder, i.e., for every realiza- 
tion of a single-particle random-matrix spectrum e^ and 
random interaction matrix elements Vap-^-yS in Eq. JJ}. 
This generates an ensemble of HF eigenvalues and HF 
eigenfunctions. We use the ensemble to calculate aver- 
age values of observables as well as fluctuation proper- 
ties (e.g., variances and/or distributions). In particular, 
we are interested in the statistical properties of the in- 
teraction matrix elements in the HF basis, and of the 
single-particle HF energies and wave functions. 

All our calculations are done for spinless electrons. 
Therefore, we use the antisymmctrized matrix elements 
defined in the previous Section. We disregard the average 
matrix elements which (for a fixed number of electrons) 
would add only a constant to the Hamiltonian. To sim- 
plify the notation we omit the superscript A. It goes 
without saying that inclusion of spin degrees of freedom 
would be highly desirable but would obviously compli- 
cate the calculations considerably. The present paper is, 
thus, an exploratory study into the use of the HF ap- 
proximation for a disordered system with interactions. 
Similar studies were previously carried out in the con- 
text of specific models, using the Coulomb interaction or 
a Hubbard-like short-range interaction^Hii^iH -^^^j^jig q^j. 
work is based on the generic many-body random ensem- 
bles introduced in Section IhI 



In Section IlII Al we discuss the HF approximation in 
the context of the induced many-body ensembles. Sec- 
tions lIH Bl and HrrCl are devoted to the study of the statis- 
tical properties of the interaction matrix elements in the 
self-consistent HF basis. In particular, in Section Fill Bl 
we derive perturbative analytic expressions for the first 
two moments of the HF interaction matrix elements (as- 
suming small u and a fixed single-particle spectrum), 
while in Section IHI CI we study numerically the statis- 
tical properties of the HF interaction matrix elements. 
Finally, in Sections IIIIDI and IIII El we study the statis- 
tical properties of the HF single-particle spectrum and 
wave functions. 



A. Hartree-Fock approximation for an interacting 
disordered system 

In the eigenbasis of ft.'-"-', the single-particle HF Hamil- 
tonian for n electrons is given by 

hai — eaSay + ^ Vap-^SPSp , (32) 

I3S 

where p is the density matrix 

n 



(33) 



Here V'z = X^a ^ii'^)\'^) ^-^^ the lowest n HF eigenstates. 
The latter satisfy 

^ hajipii^) = e["V;(") (for each a) (34) 



where e^ are the single-particle HF energies for n elec- 
trons (for n = 1 there are no interactions and e\ — £i)- 
We note that the single-particle HF Hamiltonian ha^ de- 
pends on the number of electrons n, but for simplicity of 
notation we do not indicate this dependence explicitly. 

The HF equations l|34|) are self-consistent since the HF 
Hamiltonian itself is determined in terms of its n lowest 
eigenstates ^i [see Eqs. (|5^ and (|S^ ]. The HF equations 
are solved by iteration. Starting from a particular real- 
ization of Eq and Va/s-^ys, we calculate the zeroth approxi- 
mation ho to the HF Hamiltonian, using Eq. H32() and the 
density matrix p(°) given by Eq. (123, with f/'; (/?) — Sip 
(i.e., the eigenstates of /i *■"•'). Diagonalization of ho yields 
a new set of single-particle eigenvalues ei-a- and wave 
functions V'ct (a) ■ The latter are used to construct the 
density matrix p^^^ from Eq. (|33|1 . All this yields the 
one-body Hamiltonian 



foliar — ^lia^ar 






(35) 



Diagonalizing hi we find its eigenvalues e2;cr and eigen- 
states V'ct ■ The latter are used in turn to construct a 
density matrix p(^\ etc. In general, p'^'^^ ^ p*^^-* and like- 
wise for the eigenvalues, and the iteration continues until 
the procedure converges to a self-consistent solution for 
tpi and ej" . To ensure convergence, it is often necessary 
to replace, e.g., p^^^ by p^^^ = (1 — A)p(^) + Ap^^^ where 
A (0 < A < 1) is a suitably chosen parameter. 

We define the HF basis as that basis in which the 
HF Hamiltonian is diagonal, h^r = £a S^t and par = 
Scrr'^in — a) where Q{x) = 1 for a; > and Q{x) = for 
a; < 0. The HF eigenstates a are arranged in ascending 

(n) 

order of the HF single-particle energies e^ . The latter 
can be written as 



where haJ = X^a £Q|V'cr(a)P is a diagonal matrix element 
of the single-particle Hamiltonian /i^°) in the HF basis cr. 
We emphasize that in general, ha a ^ e^- 



B. Moments of the Hartree— Fock matrix elements: 
perturbative approach 

In Section IlII CI below, we report on numerical HF 
calculations. The calculations are done for each real- 
ization of the many-particle ensemble (O , and statistics 
are collected from different realizations. This is done for 
either a fixed single-particle spectrum (where only the 
single-particle eigenfunctions are varied) or for a random 
single-particle spectrum (where both the single-particle 
energies and wave functions are varied). In either case, 
the resulting distribution of the two-body matrix ele- 
ments in the HF basis turns out to be Gaussian. We 
compare the numerical results for n fermions with the- 
oretical predictions based upon perturbation theory. To 
this end, we compute in the present Section the first and 
second moments of the two-body HF matrix elements 
perturbatively in leading order in u'^ for a fixed single- 
particle spectrum. We do so by using the statistics of 
the true (rather than of the induced) two-body ensem- 
bles introduced in Sectional This is legitimate because 
we average expressions that contain only the two-body 
matrix elements (and not other terms depending upon 
the eigenfunctions \a) oi h^°\ see Section III C|) . We con- 
fine ourselves to the case of orthogonal symmetry. We are 
particularly interested in the diagonal matrix elements. 

1. Perturbation expansion 

We recall Eqs. ^, ^ and ^. To first order in 
V, the HF Hamiltonian is obtained by substituting in 
Eq. H32|l the zeroth order expression for the one-body 

density matrix given by pgj = SspQin — [3). (We assume 
that the eigenvalues Eq. are arranged in ascending order) . 
Thus, 



ho;a-i — e-aOa-i + 2_^ '"aS-fS 



(37) 



(5=1 



where v^'^'' denotes a matrix element in the basis of eigen- 
states I a) oi h^°\ 

" (0) 

Regarding ^ '^^al-tS ^^ ^ perturbation to the Hamil- 

5=1 

tonian h^^\ we see that a HF eigenstate is (to first order 
in v^^') given by 

|a)HF = |a) + ^rW|A.), (38) 



where 



"-rr "'rrrr ^^ 



r=l 



(36) 



.(1) 



J—Y^Q[n-5)vf,.^^,. (39) 
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Using Eq. (|38() . we find that to second order in u*^"^ the matrix element I'aS-.-yS in the HF basis is given by 
I 



^^,_^(o) +\^A^)J°') +\^r^^K(") ,^Ji)A0) , V- (1) (c 






(40) 



/i^Q 



M#/9 



/^7^7 



^#-5 



To calculate the variance of a HF matrix element Va/3;7<5 to fourth order in u, we need to expand Vai3--/5 up to third 
order in v^^\ This requires an expansion of the HF state to second order in v^^\ 



\a)RF 






(41) 



where r*-^' and r"' are of first and second order in w^°\ respectively. The matrix element is then expanded as 



Vaf3-j5 



\ M#a /J5^/3 /J5!^7 P5^(5 y 



(1)^2 1 „(0) 

/3;7<5 









(42) 



where r^^ is given by Eq. H39|l . To find Tq,^ we have to expand the HF Hamiltonian to second order in v^^^ and use 
the perturbation expansion for the single-particle wave functions up to second order. The HF Hamiltonian /i'^^ is 
found to second order in v^^' by expanding the one-body density to first order in v^^\ 



Pst3 = E ^' - ^^'(^) = '^*/3®(" - /^) + np<d{n - 5) + r(,y e(n - /3)](1 - 5ps) 



1=1 



Substituting this into Eq. 1)32(1 we find 

/^L'-J = e„'5„7+E«il^©("-'5) + ^z;i°^,[rge(n-<5)+rge(n-/3)](l-5^,) 

S I3S 

= e„<5„7 + E «S7^0(" -^) + E(^S;7^ + ^'^Ls)r's'^^in - '5)e(/3 - n - 1) . 

<5 /3,<5 



For r^^) this yields 



r.(2) 



Ca - e 



— ^e(n-5) 



>? 






(43) 



(44) 



(45) 



For the diagonal matrix elements this yields 



(rW)2 I v^''^ 



I^L^a fj.^13 



2E(^i'J+^i'^>it.+«-/5 



M^a 



Z^ ' ail' Pf \"ti,y-a(3^ '^ti.p-ai 






(46) 
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and, to fourth order in w'°' 



(0) (0) 









2 V 7/"^ ,.(1)^.(1),/°) +cy^B 



-4 V v^°^ M)r^^Uv^°^ +v^°^ 

^ Z_^ "a(i;afi' a^l' pi, yJ^i^-ap^ "^l.p■av 
l_i=/,a;v^P 

° Z^ ' ati"i_ip-aP' Pv'^av;aP ■ 



^ Z^ ^aiJ.^^p-aP^aiy^„p-ap + " ^ P 



(47) 



^. yli^erage Values 

We first calculate the average HF matrix element in the ensemble (iii), the two-body Gaussian ensemble. For a 
diagonal HF matrix element Vap-ap we can rewrite Eq. H4()|l as 



VaP;aP - V^p.^p + ^ Z^ ^ailVf^p-aP + ^ 2^ ^,3/.«pa;/3a ' 



(48) 



where the third term on the r.h.s. is obtained from the second term by exchanging a and /?. To find the mean value 
of Vap-ap from Eq. H48(l . we calculate for fJ. ^ a 



''^ n5;aS''^ fj.p;ap = 2u'^{Ssp ~ &bp^aP^ba + & iLp^&ii^aP^&a ~ 6fj,p65fj,6sp 

+^^tQ^(5/3 — S^pSsa^afi + (>i_ip6saSapSsp, — S^a^SpSapSsfj.) 

= 2u {Ssp - SspSap - Sf_ipSsf_i6sp) = 2m Ssp{l - Sap - Sf_ip) , 
where all other terms vanish since /x 7^ a. Averaging Vap-ap in Eq- I48|) using Eq. (|49|l . we find 



VaP;aP = V^p-aP + ^^^ C^ - ^^p) 



0("-/3)(Er^ 



vM#a 



Eq — £^l £« — £,3 



+ a^ f3 



4u2(l - ,5„^) 



e(n-a) J2 -^— + e{n-P) ^ — 



/./a,/3 



ep~e 



li^a.p 



^a t/i 



For a general matrix element t'a/3;75 we find 



(49) 



(50) 



Vap-jS ^ ^ ralvfp.^^ + (a ^ /3, 7 ^ (5) + (a ^ 7, /3 ^ 5) + (q ^ 5, /? ^ 7) 
where we have used Eq. (|4()() and the symmetries of v. We have 



E-^^^S.. - -^Y.^{n-a)v^l 



(0) „(o) 



ii^a 



,(0) „(o) 



Using Eqs. lO and (HU, we find 



•fia-aaV^p.^g = 2u 5cr/3 (1 - (5^/3) ('^a7'^o-5 - '^QS'Jcr^) 



VaP;-iS = {Sa-i^pS — 5a&5 pj)Vap-aP , 



(51) 



(52) 
(53) 



where Vap-ap is given by Eq. H5U|) . 
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A similar calculation for the orthogonal ensemble (ii) yields for the average HF matrix element 



fQ/3;Q/3 = 6u^(l - Sap) 



e(n-a) J2 -^— + e{n~f3) ^ — 



li^aJ3 



e/3 



^5^a,/3 



^a ^Li 



(54) 



Comparing this with Eq. (|50l) we see that the result is the same as for the two-body Gaussian ensemble (iii) except 
for an overall factor of 3/2. 



3. Variance 



The expectation values of each of the seven terms on the r.h.s. of Eq. (|47|l can be calculated using Wick's theorem 
and Eq. (j^ . The calculation proceeds along lines similar to those of the previous Section. We do not give any details 
of this rather lengthy calculation and confine ourselves to the result. For the variance of the diagonal HF matrix 
elements we find 



<j'^iVaf3;al3) = 4m^(1 - Sa^) 



l^iu'Q{n-a) Y, I ^ -^u'Qin-P) ^ ^ 



^L^a.,|3 



(e/3 - e^) 



Ai#Q,/3 



(ea - e^)^ 



r 



(55) 



We observe that in contrast to the expressions derived 
for the average values of the matrix elements, Eq. H55(l 
cannot be averaged over the GOE spectrum. This is 
because the integrals over the eigenvalues possess a log- 
arithmic singularity. We shall return to this point when 
we compare the numerical results with the perturbative 
expressions in Section llllCBI 



Statistics of the Hartree-Fock matrix elements: 
numerical results 



We now report on numerical results concerning the 
statistics of the interaction matrix elements in the self- 
consistent HF basis, using the true random-matrix en- 
sembles of Section ^ We focus attention on the diago- 
nal matrix elements v„r;aT- Unless otherwise stated, we 
use in this Section the orthogonal "non-local" two-body 
ensemble (I22|l with m = 30 single-particle states and 
n = IQ electrons. 

In the simulations, we construct the ensemble (|2J as 
follows. We generate an m x tti GOE matrix h^^^ and 
diagonalize it to find its spectrum ea- We also gener- 
ate Gaussian variables Vap-^-yS that are uncorrelated with 
the e^'s and satisfy Eq. j^UJ). We then compute the cor- 
responding antisymmetrized interaction matrix elements 

from V^p.^s = '"af3-nS - Val3;5f 

The parameters of the ensemble are A and u. We have 
chosen the variance of the off-diagonal matrix elements 
of the GOE matrix /i(°) to be cr'^{h^°p) = Sm/An'^, so that 
the mean level spacing in the middle of the spectrum is 
fixed at A = VS/S « 0.866. For each of the quanti- 
ties presented in the figures of this Section we have used 
10, 000 realizations. 

In some of the calculations presented below, the single- 



particle spectrum is kept fixed and an ensemble of 
antisymmetrized two-body interaction matrix elements 
v^a-js is generated as discussed above. 

We used a rather low-dimensional GOE matrix to 
generate the one-body ensemble. An alternative choice 
would have been to use the center section of a high- 
dimensional GOE matrix as our starting point. Our 
choice has the obvious disadvantage of being not free of 
edge effects. However, the comparison made below of 
some of our results with those obtained from a picket- 
fence spectrum (constant spacings) shows that edge ef- 
fects are small for the quantities we investigate. Our 
choice has the obvious advantage of keeping the dimen- 
sion of all relevant matrices small. This is very useful in 
view of the number of realizations we have chosen, and 
of the need to perform a HF iteration for each of these 
realizations. 



1. Distributions 

We first consider the case of a fixed single-particle 
spectrum, allowing only the single-particle wave func- 
tions and, thus, the two-body matrix elements to vary 
randomly. For each realization of the two-body matrix 
elements, we solve the HF equations and compute the di- 
agonal matrix elements Vap-ap in the HF single-particle 
basis. 

Fig- d shows the distributions of three diagonal matrix 
elements Vap-^p for n = 10 electrons and u — 0.1. For 
a — 5,P = 5 both levels are filled, for a = 5, /? = 20 one 
level is filled, the other is empty, and for a — 20, (3 = 20, 
both levels are empty. The distributions are all well de- 
scribed by Gaussians (solid lines). The same holds true 
for the distributions generated by choosing both single- 
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A a = 5, /? = 5 
□ a = 5, /? = 20 
o a = 20, /? == 20 




0.75 



^Q/3;a/3 



FIG. 1: The distributions P(vai3-ai3) of three diagonal matrix 
elements in the HF basis for a fixed single-particle GOE spec- 
trum in the ensemble 12211 with u = 0.1. There are n = 10 
electrons. The symbols denote the results of the numerical 
simulations. The solid lines are Gaussian fits. 



particle energies and two-body matrix elements at ran- 
dom. It is, thus, sufficient in all cases to consider only 
the first two moments of Vafi^ap- 



2. Average HF matrix elements 

The top two panels (a) and (b) of Fig. [21 show the av- 
erage values of some diagonal matrix elements calculated 
for a fixed single-particle spectrum. This fixed spectrum 
was chosen as a picket-fence spectrum (equal spacings) 
in panel (a) of Fig. |2 and as a GOE spectrum in panel 
(b) of Fig. 121 The perturbative expression (|50|l shown by 
solid, dashed and dashed-dotted lines is in good agree- 
ment with the numerical results. The agreement is better 
for a picket-fence spectrum than for a GOE spectrum. 

We turn to the general case where the statistics are 
collected by varying both the single-particle GOE spec- 
trum and the single-particle wave functions. Such fully 
averaged HF matrix elements are shown in panel (c) of 
Fig.El 

According to the perturbative expression (|50|l , the full 
ensemble average of a HF diagonal matrix element is 
given by 4u^(l — 5ai3)Cap with 



Cap 



Q{n~a){ J2 



\t,^a,f3 



e/3 - e^ , 



e(n-/5)( ^ 



\ai#q,/3 



(56) 



The symbol (. . .) denotes the GOE average. In Fig.|21we 
show a three-dimensional plot of Cap versus a and /3 for 
m = 30 and n = 10. 




FIG. 2: Average matrix elements Vap-ap in the HF basis ver- 
sus a for fixed /3. We show results for three values of /3 
(/3 = 6, 12, 25). We use the ensemble ^ with m = 30 and 
for n = 10 electrons, (a) The data are generated for u = 0.1 
and for a fixed picket-fence single-particle spectrum (equal 
spacings). The symbols give the numerical simulations, and 
the lines give the perturbative result 1501 . (b) As in panel 
(a) but for u = 0.05 and for a fixed single-particle GOE spec- 
trum, (c) The results are for the same ensemble as in (a) , but 
the average is taken over the full ensemble 12211 . i.e., over both 
the single-particle GOE spectrum and the two-body interac- 
tion matrix elements of Eq. Q. The symbols represent the 
numerical simulations and the lines the perturbative expres- 
sion 1501 . averaged over the single-particle GOE spectrum. 



For comparison, we show in Fig. ^ the full ensemble 
average of a diagonal HF matrix element (averaged nu- 
merically over both the two-body matrix elements and 
the single-particle GOE spectrum) versus a and f3 for 
u ~ 0.1. The result agrees well with Gu^Cap if we take 
Cap from Fig.El 



3. Variance 

The two-body matrix elements which serve as input 
in the HF calculation have a Gaussian distribution with 
second moments defined in Section IIIB II Our numeri- 
cal simulations indicate that the variances a^ of the HF 
matrix elements are very close to these input values. We 
observe that in the HF basis a decreases slightly, the 
size of the decrease depending on whether the diagonal 
matrix element corresponds to fiUed-fiUed, fiUed-empty 
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FIG. 3: Three-dimensional plot of Ca/s [see Eq. I|56|l ] versus a 
and /3. The results shown are for n — 10 electrons occupying 
m — 30 orbitals. 




filled-filled matrix elements. 

We have also compared our numerical results for a fixed 
choice of the single-particle spectrum with the perturba- 
tive expression (|55l) . This expression seems to fail no 
matter how small u was chosen. (When u becomes too 
small, the statistical noise prohibits a meaningful com- 
parison). While the numerical result is always close to 
4w^, this is not true for the expression (|55|l which at times 
may even become negative. We speculate that this fact 
may signal a breakdown of perturbation theory but have 
not followed up this point. 



b 
I 




0.05 



0.25 



FIG. 5: Top panel: The standard deviation of a HF matrix 
element a{vai3-ai3) versus u for the "local" interaction ensem- 
ble 11171 . We show results for three different types of diagonal 
matrix elements, corresponding to filled-filled, fiUed-empty 
and empty-empty single-particle levels. The results shown 
for each type are averaged over all matrix elements of that 
type. Bottom panel: The fractional deviation 1 — a^ /Gu^ of 
the variance a^ of the matrix elements in the HF basis from 
their variance 6m^ in the non-interacting basis versus u. Same 
conventions as in the top panel. 



FIG. 4: Three-dimensional plot of the full ensemble average 
of Vai3-ai3 vcrsus Q and /3. We use the "local" interaction en- 
semble with u = 0.1 for n = 10 electrons occupying m — 30 
orbitals. The results shown here agree well with the pertur- 
bative expression Qv?Cai3 where Ca/s is shown in Fig.|3| 



or empty-empty single-particle orbitals. There are fluc- 
tuations because of the finite size of the ensemble. To 
improve the statistics we further average the variances 
separately for these three different types of matrix ele- 
ments. The top panel of Fig.|31shows the resulting stan- 
dard deviation <7{vaf3-ai3) for the "local" interaction en- 
semble (|17|l versus u. The input value is a = vGu. The 
difference between the HF results and the input value is 
amplified in the bottom panel of Fig. \^ The variance 
of the empty-empty matrix elements is practically un- 
changed, while the largest decrease is observed for the 



D. Statistics of the Hartree— Fock single— particle 
energies 

The statistical results presented in the following are 
calculated for the entire ensemble, i.e., by collecting 
statistics over both the single-particle spectrum and the 
single-particle wave functions (i.e., the interaction ma- 
trix elements). Unless stated otherwise, the number of 
realizations used in the numerical simulations is 10, 000 
(for each value of u) . 

The mean level spacing A of the single-particle GOE 
levels is given approximately by the inverse semi-circle 
law and increases towards the edges of the spectrum. 
We now consider the mean level spacing Ahf of the HF 
single-particle spectrum calculated by averaging over the 
entire ensemble. Not surprisingly, we find that Ahf in- 
creases monotonically with u. Fig. shows Ajjf versus 
the level index a for different values of u. The solid line 
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is the GOE result. We notice the large peak at n = 10. 
It is due to a large spacing between the highest filled and 
the lowest empty HF level. This is a well-known feature 
of HF calculations and is referred to as the HF gap. We 
notice also the slight increase in the mean spacing close to 
the HF gap. In our discussion of the statistics of the HF 
single-particle energies Ca , we treat filled levels, empty 
levels and the HF gap separately. 
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FIG. 6: The mean level spacing Ahp of the single-particle HF 
levels versus level index a for u = 0.05 (triangles), u = 0.1 
(squares) and u = 0.15 (circles). The results shown are for 
n = 10 electrons and m = 30 single-particle states. The 
solid line describes the mean level spacing A for a GOE of 
dimension m — 30. Inset: the mean level spacing in the 
center of the spectrum versus u. 



1. Filled levels and empty levels 

Let s denote the spacing of a pair of successive filled or 
empty levels, measured in units of the local mean level 
spacing. The nearest-neighbor spacing distribution P{s) 
is found to follow the GOE Wigner distribution. We 
observe such behavior even for the pair of filled (empty) 
levels closest to the HF gap. This is demonstrated in 
Fig. [3 where the numerical distributions (histograms) for 
P{s) are compared with the Wigner distribution (dashed 
lines) for different values of u and for both the highest 
two filled HF levels (top panels) and the lowest two empty 
HF levels (bottom panels). 



2. Hartree-Fock gap 

We simplify the notation by writing Vap = Vap-^ap- In 
the HF basis and for n electrons, the HF gap is given by 




FIG. 7: Nearest-neighbor spacing distribution P{s) versus 
spacing s where s is measured in units of the local mean level 
spacing of the corresponding HF levels. Top panel: The his- 
tograms show P{s) for the highest two filled single-particle 



HF levels e 



(n) 



and e„ (with s — (e„ - 



,(") 



i)/(' 



_(n) 



-i"-'l)), 



and for three values of u: (a) u = 0.05; (b) u = 0.1; and 
(c) u = 0.15. The dashed lines give the GOE Wigner dis- 
tribution. Bottom panels: same as in the top panels but for 
the two lowest empty single-particle HF levels ejj'^i and e„'^2- 
The values of u are (d) u — 0.05; (e) u = 0.1; and (f) u = 0.15. 
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-h)al)+Vn+l,n+}XVn+l,T-Vn,T) 



= 1 



(57) 

The squares in the top panel of Fig. |Slshow the numer- 
ical results for the average gap, calculated from the full 
ensemble © with a "local" interaction [Eq. (|17|l ]. To es- 
timate the average gap, we use in Eq. l(57|l the perturba- 
tive expression H54|l for the diagonal HF matrix elements 
which we average over a single-particle GOE spectrum 
(see Caj3 in Section IIII CI) . We also need to know the 



,(0) 



,(0) 



which arc found to be 



average values of h^^i „_|_i 
0.895,0.868,0.836, and 0.765 for u = 0.05,0.1,0.15 and 
0.25, respectively (the corresponding value in the absence 
of the two-body interaction is A = 0.908). The result- 
ing estimate for the average gap depends quadratically 
on u, and is shown by the dotted-dashed line in the top 
panel of Fig. |S1 For small values of u (< 0.15) our es- 
timate provides a good approximation for the average 
HF gap, but for larger values there are deviations, in- 
dicating the breakdown of the lowest-order perturbative 
approach. We note that the average interaction Q con- 
tributes a constant vq to the average gap. In our calcu- 
lations we have set the average interaction to zero, and 
this constant (but large) contribution to the average gap 
is not shown in our numerical results. 

The distribution P{s) of the HF gap, shifted by its av- 
erage value, is shown in Fig.|51for several values of u (his- 
tograms) . In Section IIV Al we motivate an approxima- 
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FIG. 8: Top panel: average HF gap versus u. The squares are 
the results of simulations of the full ensemble Q. The dotted- 
dashed line is obtained by using Eq. 15711 and a spectral GOE 
average of the perturbative expression H5()|l . Bottom panel: 
the squares depict the width cr of a Gaussian the convolution 
of which with a Wigner distribution is fitted to the gap dis- 
tributions in Fig. |5| The dotted-dashed line is \/6w, and the 
triangles (connected by the dashed line) are the standard de- 
viations (j{vn+i,n) of the corresponding HF interaction matrix 
element (see text). 



tion to this gap distribution as a convolution of a (shifted) 
Wigner distribution with an appropriate A, and a Gaus- 
sian. The dashed lines in Fig. |51 are the corresponding 
Wigner distributions, while the solid lines describe the 
convolutions for which the width a of the Gaussian is fit- 
ted. This fitted value of a is shown in the bottom panel 
of Fig. 13 versus u (squares). The dotted-dashed fine in 
Fig. |S1 is the input value \/&u for the standard deviation 
of a diagonal matrix element. The triangles give the stan- 
dard deviation a{vn+i,n) of the diagonal matrix element 
Vn+i^n in the HF basis for n electrons. We observe that 
a is somewhat smaller than the standard deviation of the 
corresponding HF interaction matrix element. 



3. Spacing correlator 



In Section IIIID II we have seen that the nearest- 
neighbor spacing distributions of both, the occupied and 
the empty HF states, follow RMT statistics. But what 
about the correlation between occupied and empty lev- 
els? In this Section we show that the HF gap modifies 
that correlation so that it differs from the RMT predic- 
tion. To that end, we define the spacing correlator 



An) 



SaSa+2 — SaSa+2 

a{sa)cr{sa+2) 

Jn) 



(58) 



where s„ = ea^' — e^^l^ (a = 2, ...,7n) are nearest- 
neighbor spacings and cr(sQ,) is the standard deviation 
of Sq. For the HF ensemble, the correlator Cr. measures 



FIG. 9: The calculated distribution of the shifted HF gap 
(histograms) for four values of u: (a) u — 0.05; (b) u — 0.1; 
(c) u = 0.15; and (d) u — 0.25. The gap is shifted to have the 
average value zero. The solid lines describe convolutions of a 
shifted Wigner distribution (with an appropriate A) with a 
Gaussian. The width of the Gaussian is fitted. The dashed 
lines are the shifted Wigner distributions. 



the spacing correlation across the gap, i.e., the correla- 
tion between the highest spacing of the filled levels and 
the lowest spacing of the empty levels. 

For the GOE, the correlator l|58|l is denoted by C{r = 
l;t = 0) and explicitly given in Ref. 38. The correla- 
tor is negative and approximated by C'(r = l,t = 0) ~ 
— l/[47r^(4/7r— 1)] « —0.093. A more accurate numerical 
value is -0.087. 

In the HF approach, the correlator (|58|) is small and 
exhibits large sample-to-sample fluctuations. Thus, a 
large number of realizations is required to obtain a re- 
liable value. Table I shows c„ of the HF single-particle 
spectrum for several values of u. The results were ob- 
tained using 50, 000 realizations (for each u) of the "local" 
interaction ensemble with m = 30. We observe that the 
correlator is negative and its magnitude becomes smaller 
as u increases. Hence the HF gap has the effect of weak- 
ening the correlation between filled and empty levels. 



u 


correlator c„ 





-0.084 ±0.011 


0.05 


-0.078 ±0.011 


0.10 


-0.061 ±0.012 


0.15 


-0.048 ±0.011 


0.25 


-0.027 ±0.013 



TABLE I: The spacing correlator c„ across the HF gap for 
several values of u (using 50, 000 realizations for each u) . 

On the other hand, when both spacings are between 
filled states (a < n — 2) or both between empty states 
(a > n± 2), we find that c^ is close to its standard GOE 
value (except for edge effects). This provides additional 
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support for the hypothesis that the filled and the empty 
HF states separately satisfy RMT statistics. 

Fig, irni demonstrates the behavior of the correlator Ca 
as a function of a for u = 0.1 and u = 0.25. For a < n — 2 
and a > n + 2, Ca fluctuates around its RMT value. 
However, the correlator c„ between spacings on opposite 
sides of the gap decreases in magnitude with u, while the 
correlators c„_i and c„+i (for which one of the spacings 
is the gap itself) are enhanced in magnitude compared 
with the RMT value. 



0.1 

0.05 



-0.05 

-0.1 

-0.15 



by 



f-fli- 




-i4-J-t-rii 



5 10 15 20 25 

a 



FIG. 10: The spacing correlator Ca [see Eq. 1)58^ ] versus a for 
the HF spectrum of the "local" interaction ensemble with u = 
0.25 and tti = 30 (values of a near the edges of the spectrum 
are omitted). The results shown are calculated for 50,000 
realizations of the ensemble and include statistical errors. The 
dashed line is the GOE value of the correlator —0.0863 ± 
0.0008 calculated from 500,000 realizations and averaged over 
a. Inset: the correlator Cn across the gap versus u. 



E. Statistics of the Hartree-Fock single— particle 
wave functions 



The induced ensemble is invariant under orthogonal 
(unitary) transformations of the single-particle basis. We 
show now that the corresponding single-particle HF en- 
semble is also invariant under orthogonal (unitary) trans- 
formations. For simplicity, we present the proof for the 
orthogonal case only. 

We consider a particular realization (0) of the in- 
duced ensemble in a fixed basis \i). Under an orthog- 
onal transformation O of the single-particle space with 
a] —> OajO = ^j, a[,Oiii, we obtain a new realization of 
the many-body ensemble. 



H^H = OHO-' = ^ hf^a\a, + - ^ itrM^^^] 



ij;kl 



aiak 



(59) 



where the new single-particle Hamiltonian h^^' and in- 
teraction matrix elements v are related to w-^' and v 
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''ik 






Y.O,,hf{0-')ik- 



]l 



Y, 0.,eO,,,v^^,.^,,i,iO~%,k{0~'),i .(60) 

i' j' k'l' 



The Hamiltonians ft-'-*'-' and h^^^ have identical single- 
particle spectra but different single-particle eigenstates. 
Both are members of the single-particle GOE. Likewise, 
the Hamiltonians H and H are both members of the in- 
duced many-body ensemble. 

In the fixed basis the single-particle HF Hamilto- 
nian (|32(l has the form 



riik - ri-^ 






ij;kpPpj 



(61) 



Here ppj = {j\p\p) = J2i^'iip)^iiJ) is the matrix rep- 
resentation of the single-particle density operator p = 
ELi \i^i){M- As in Eq. §^, IV'/) denotes the l^'' HF 
single-particle eigenstate, ipi{j) its projection onto the 
fixed basis state j, and the sum is over the lowest n 
HF levels. We consider the transformed single-particle 
Hamiltonian 



hik = YOi-jhji{0 '] 

ji 



Ik 



(62) 



To prove the orthogonal invariance of the ensemble of HF 
single-particle Hamiltonians, we have to show that h is 
the HF Hamiltonian for the new realization H59(l of the 
many-body ensemble. Using Eq. H61|l we can write h as 



hik = kl +J2^ti;kiP^i 



(63) 



where 



Opo-' = J2 \^i){i'i\ : with \i,i) = d\^i[ 
1=1 



(64) 



This shows that hik is indeed the HF Hamiltonian of the 
realization (|59|l provided that p is the density matrix of 
the lowest n eigenstates of h. But Eq. 1)62(1 shows that 

(n) 

h has the same eigenvalues e^ as h and that the eigen- 
vectors of h are given by the second of Eqs. 1(64(1 . This 
completes the proof. We conclude that the invariance 
of the many-body ensemble under orthogonal transfor- 
mations implies the orthogonal invariance of the HF en- 
semble. In the GOE, the statistics of the eigenvectors are 
determined by the orthogonal invariance of the ensemble. 
Thus we expect the eigenvectors of the HF ensemble to 
obey the same statistics. For example, the distribution 
of the components V'/(i) of the HF eigenvector |'0/) in 
the fixed basis of states \j) must depend on orthogonal 
invariants only. For a single eigenvector the only such 
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orthogonal invariant is ^ i^fij)- Since the eigenvector 
is norniahzed to unity, the probabihty density is 



PiM^),M2),---)^s{J2^fij) 



(65) 



As for the GOE, integrating over aU components except 
the j-th, we find 



p{Mj)) = ^' 



-1/2 



r(2 



r(^; 



[i-V^f(j)]^ 



(66) 



where m is the dimension of the single-particle space. In 
the limit of large to, P{i^i{j)) can be approximated by a 
Gaussian, P{Mj)) « e-'"'^?(^')/2. 

We have computed numerically the statistics of an HF 
eigenfunction component in the "local" interaction or- 
thogonal ensemble (|17|l . Since there is an ambiguity of 
the overall sign of the wave function, we have calculated 
the distribution of y = 'rmpf{j) (the normalization is cho- 
sen to satisfy y — 1). It follows from H66|) that 



P{y) = vr 



-1/2 



r(^) 



r(^) 



(-(TO,)-V2 i_A ^ 



(67) 



In the limit of large m, this is just the Porter-Thomas 
distribution 



P{y) = (27ry) 



-i/2g-y/2 



(68) 



Fig. 111! shows numerical results for the distribution 
of P(lnj/) versus Iny for an HF eigenstate in the middle 
of the spectrum and for m = 30. Since the distribution 
is independent of the particular component, we collect 
statistics from all components. Results are shown for 
u = 0.05,0.1,0.15 and 0.25 (symbols). For all values of 
u the distributions are well described by Eq. (|57|l (solid 
line). For reference we also show the Porter-Thomas 
distribution (dashed line). The observed small deviation 
from the limiting case of a Porter-Thomas distribution 
is a finite-size effect. 

As in the GOE, it is also possible to calculate from 
Eq. (|65|) the correlator for two different squared compo- 
nents of an HF eigenfunction. 






1 



(69) 



where a'^{tpf{i)) is the variance of ipfii)- This correlator 
is rather small and in the numerical simulations exhibits 
large fluctuations. However, on average it agrees with 
the analytical result Eq. (|69|l . 

In Sections IIII Dl and IIIIEI we have thus shown that 
the fluctuation properties of the HF ensemble are very 
close to those of RMT: The single-particle HF wave func- 
tions obey the same statistics, and the single-particle HF 
eigenvalues do, too, except near the HF gap which plays 
a special role in the spectrum. The correlations between 
filled and empty levels are weakened by the gap. 




FIG. 11: The distributions P{lny) versus Iny, where y is 
the square of an HF wave function component (scaled to give 
y — 1) for u — 0.05 (triangles), u — 0.1 (squares), m = 0.15 
(circles) and u = 0.25 (diamonds). We use the "local" inter- 
action orthogonal ensemble 1171 with m — 30 and n = 10. 
The solid line is the finite-m RMT distribution 16711 and the 
dashed line is the Porter-Thomas distribution 16811 . A section 
of the graph is magnified to render more details. 



IV. HARTREE-FOCK APPROACH TO 

DISORDERED SYSTEMS: ADDITION OF 

ELECTRONS 

In the HF approximation, the many-particle ground- 
state energy -Ehf ("■) of n electrons can be written in sev- 
eral ways. In an arbitrary fixed single-particle basis, we 
have 

EnFin) = ^ h^ajPia + -^ X! PiaVa(3-^SP5l3 , (70) 

where p is the self-consistent density matrix of Eq. H33|l. 
An alternative expression (in terms of the single- 
particle HF energies) is 



n 



Val3-j5P5[3 



(71) 



fjS 



where the double counting of the interaction terms (in 
the sum of the single-particle HF energies) is corrected 
by the second term on the r.h.s. of Eq. H71(l . Both of 
these expressions simplify when written in the HF basis 
[in which p^a — Sa-yQiji — a)] 

n -. n 

a=l a.J3=l 

n -J n 

= E^i"^-2 E '^"ft"^- (72) 

a=l a,/3=l 
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A. Peak spacing distribution 

The position (as given by the gate vohage) of a 
Coulomb-blockade peak at low temperature is indicative 
of the change in the ground-state energy of a quantum 
dot due to the addition of an electron (known as the 
addition energy). The spacing A2 between successive 
Coulomb-blockade peaks (for short: the peak spacing) is 
then given by the second-order difference of the ground- 
state energies versus particle number. In the HF approx- 
imation we have 



A2 ~ EiiF{n + 1) + ^HF(n - 1) - 2^HF(n) 



(73) 



In the CI model, i.e., when only the average interaction 
is taken into account, A2 = (cn+i — Cn) + vq, and the 
distribution of A2 is a shifted Wigner distribution. Ex- 
perimentally, the distribution is closer to a Gaussian(S*S*i 
and this was understood to be a residual interaction ef- 
fect. 

The histograms in Fig.[T21show P(A2 - A2) in the HF 
approximation for the same cases as in Fig. ^ The solid 
lines are again convolutions of a Wigner distribution with 
a Gaussian of width a for the spacing distribution (this 
will be motivated in the following Section). The fitted 
values of a are shown versus u as squares in Fig. 1131 The 
dotted-dashed line, triangles, and dashed line are the 
same as in Fig. [S] We note that the average interaction 
contributes to A2. 

In the spinless case, the average interaction is given by 
Eq. |(7J), contributing a constant vq to A2, which is can- 
celed out in A2 — A2. Thus the distribution P(A2 — A2) 
is not affected by the average interaction. This would no 
longer be true if spin were included: The average interac- 
tion in Eq. (j^ has additional terms (i.e., an exchange in- 
teraction, and, in the orthogonal case, a Cooper channel 
interaction), which are not constant for a fixed number 
of electrons. 




0.25 



FIG. 13; The squares depict the width cr of a Gaussian the 
convolution of which with a Wigner distribution is fitted to 
the peak spacing distribution P(A2 — A2) in Fig. 1121 Dotted- 
dashed line, triangles, and dashed line are as in Fig. |H| 



B. Koopmans' approach 

As an electron is added to the quantum dot, the 
self-consistent single-particle HF wave functions are ex- 
pected to change. In Koopmans' limit, this change is 
neglected^ This assumption implies predictions for the 
peak-spacing distribution which we now test. 

We first calculate the addition energy £'hf("- + 1) ~ 
E-tiF{n) in Koopmans' limit. We use the HF wave func- 
tions for n electrons to write expressions for E'hf (n) and 
-Ehf(^i + 1). For £'hf(?i) we use the exact HF relation 
Enpin) = Y^a^ihaa + ^Y.lf3=i^ai3- Usiug the same 
single-particle wave functions, we write for EnF^n + 1) 
the approximate relation Enpin -I- 1) « ^^ 
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n+l 



,n+l , (0) 
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We obtain 
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FIG. 12: Peak spacing distribution P(A2 — A2) for the same 
set of u-values as in Fig. |5] 
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(74) 
A similar expression can be obtained for the n—l—i-n 

transition, i.e., _Ehf("-) ~ -£'hf('t- — 1) ~ £« , where we 
have used the HF wave functions of n — 1 electrons. We 
then find from Eq. (|73|) the approximate relation^ 



A, 



-n+l 



^(«-l) - (J^^ 



.("-1)^ 



\'^n+l ''n+l / ^ \'^n+l '^n ) 



Vn+l,n + S , 



(75) 



where s = e. 



in-l) _Jn-l) 
n+l 



is the spacing between the two 



lowest empty levels for n — 1 electrons. In deriving the 
second approximate relation in Eq. (TSJ, we have also 

used e"_li « ^n+i ) + Vn+i,m an expression that follows 
from Koopmans' limit for the n — 1 —>■ n transition. 

An expression analogous to Eq. ifT^ can be derived 
from the n -\- 1 —^ n transition, assuming the single- 



20 



particle HF wave functions do no change upon the re- 
moval of an electron from the dot. Using the HF wave 
functions of n + 1 electrons, we find 



EnFin + 1) - EuF^n) « e 



(n+l) 
n+1 



(76) 



Applying a similar relation for the n —f n — 1 transition, 
we find another expression for A2, 



A. 



(n+l) 



.(") 



S + Vn+l^,i 



Jn+l) 



) + (4' 



(n+l) 



,(«) 



) 

(77) 



where s — e 



(n+l) 
n+l 



e„ is the spacing between the two 
highest filled levels for n + 1 electrons. 

If the HF levels and HF interaction matrix elements 
are uncorrected, Eqs. lf75|) and lf77|) suggest that the dis- 
tribution of A2 can be described by a convolution of a 
Wigner distribution (with A given by the average spac- 
ing of the two lowest empty or two highest filled levels 
in the dot, respectively), and a Gaussian whose width 
a is given by the standard deviation of the HF matrix 
element t'n+i.n- The results shown in Fig. 1121 confirm 
that the peak spacing distribution is well described by a 
convolution of a Wigner distribution with a Gaussian (A 
is determined from a mean level density that excludes 
the gap). Furthermore, for small values of u (< 0.15), 
the width of the Gaussian is well described by a{vn+i,n) 
(see Fig. I13|l . For larger values of u, we observe some 
deviations. 

In view of Eqs. (|75|l and (|77|) . it is interesting to com- 
pare directly the peak spacing distribution P(A2) with 

the distributions P(ei++/^ - ei"^) and P(ei"^^i - ei"^'^) 
using the exact self-consistent single-particle HF ener- 
gies (i.e., not using Koopmans' limit) ii^ The results are 

shown in Fig. [H for u = 0.1. Both P(ei++^^ - e^"^) and 



P{' 



-n+l 



.("-1)^ 



approximate f (A2) very well. 



If we use Koopmans' limit simultaneously for the addi- 
tion and removal of an electron in an n-electron dot, i.e., 
for the transitions n ~f n + 1 and n — > n — 1, we obtain 
the relation A2 ~ e„_li — e„ . The quantity e^_li — £« 
is just the gap in the n-electron dot. This fact motivates 
the description of the gap distribution in Fig.|51as a con- 
volution of a Wigner distribution with a Gaussian. The 
width (T of the Gaussian (found from the gap distribu- 
tion) is rather close to cr(wn+i,n) (see Fig. ISJ, although 
not as close as the value of a found from the peak-spacing 
distribution (see Fig. lT^ . 

In Fig. 1141 we also compare the gap distribution 



Jn) 



An) 



P{^n+i ~ ^n ) with the peak-spacing distribution (with- 
out subtracting the average values of the respective quan- 
tities). We observe that the gap distribution is similar in 
shape to P(A2) (both are convolutions of a Wigner dis- 
tribution with a Gaussian) but is shifted to the right. 

This suggests that the approximation A2 ~ ^n+i ~ ^n 
does not work as well as Eqs. (|75|) and (|77|l . in particular 
for the average values. 



The qualitative difference in the various approxima- 
tions for A2 can be explained as follows. The single- 



particle HF energies e„"^i and en in Eq. (|75|) corre- 
spond to empty levels and thus contain both HF matrix 
elements connecting empty and filled orbitals. The aver- 
age values of these HF matrix elements are similar (see, 
e.g., in Fig. Inland Eq. (|54|l ) and cancel out when their 
differences are taken to find A2. Similarly, the HF ener- 
gies e^+i and £« in Eq. H77|l both correspond to filled 
levels. Both contain HF matrix elements between two 
filled levels whose averages are again similar and cancel 
out upon taking their difference. Because of this can- 
cellation effect, the approximations in Eqs. H75|l and H77|) 
work better than those in Eqs. H74|) and (|76|l . respectively. 
However, the gap is the difference between an empty level 
and a filled level {e„J_i and e„ , respectively). The HF 
interaction matrix elements are of the type empty-filled 
for the empty level and filled-filled for the filled level. 
The average values of these different types of matrix ele- 
ments are rather different and their differences contribute 
to the large average HF gap observed in Figs. ISl and 1141 
We note that both A2 and the gap include a large addi- 
tional constant vq (mostly charging energy) which is not 
shown in our numerical results (since we set the average 
interaction in the original ensemble to zero). 



Qh 0.4 




FIG. 14: The peak spacing distribution -P(A2) (left his- 
togram) is compared with P(e^"3L 



A")\ 



(dashed line), 

P{e„^i — e„~ ) (dotted-dashed line) and the gap distribu- 
tion (right histogram). The ensemble has u — 0.1. 

Another motivation for comparing the peak-spacing 
distribution with the above three distributions (using in 
each case the exact HF levels for the appropriate number 
of electrons) is provided by inequalities for the addition 
energy that are exact in the HF approximation. The 
basic inequalities arei^ 



(n+l) 
^n+1 



<EiiF{n + l)-EiiF{n)<e 



(n) 



(78) 



Both inequalities follow from the variational principle for 
the HF ground-state energy. To obtain the right inequal- 
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ity, we use the variational principle for n + 1 electrons. 
As a trial wave function we choose the Slater determi- 
nant of the lowest n + 1 HF wave functions found in the 
exact HF solution for n electrons to obtain Eiip{n+ 1) < 

X]a=i ^aa + \ X]"/3=i ^a/3- Combined with the exact re- 
lation £;hf("-) = YTa=i hal + \ YJap=i ""^p, wc find the 
right inequality in H78|) . Similarly, the left inequality in 
(|78|l is obtained by using a trial Slater determinant com- 
posed of the lowest n orbitals found in the exact HF 
solution for n -I- 1 electrons to put an upper bound on 

£'hf('^)- 

Rewriting H78|l with n replaced by n— 1, we obtain HF 
inequalities for the "removal" energy _Ehf ("■—!) ^-E'hf('t-) 




-ei""'^ < ^HF(n - 1) - ^hfH < 



,(") 



(79) 



The peak spacing A2 is obtained by summing the addi- 
tion and removal energies in Eqs. (|78|l and (|79|) . We find 
that the gap is an upper bound for A2, i.e., 



A2 <e 



n+l 



.(") 



(80) 



This relation is consistent with the shift to the right of the 
gap distribution relative to the peak-spacing distribution 
(see Fig. [n|) . However, inspecting the inequalities in lf75|) 

and ifT^ . we conclude that the quantities e„' 



+1 



An-l) 



and t^j^i — en are neither upper nor lower bounds for 
A2. As discussed above, they turn out to be rather good 
approximations for A2 (for u < 0.15). 



The gap and peak— spacing distributions in the 
small u limit 



In Sections IIII D 21 and IIV Al we have shown that the 
gap and peak-spacing distributions can be approximated 
as convolutions of a Wigner distribution (evaluated with 
the HF mean-level spacing) and a Gaussian whose width 
a is close to cr(w„+i^„) (for the "local" interaction ensem- 
ble a{vn+i,n) ~ v6w). On might have chosen another 
approach: One may attempt to use Eq. H57() directly to 
infer the distribution of, e.g., the gap in the limit u — > 0. 
The various diagonal HF interaction matrix elements in 
Eq. (|57|l are approximately uncorrelatcd Gaussian vari- 
ables with variance cr^(fn+i,n) ~ 6u^ each. Furthermore, 



to first order in the interaction, we have ha.a ~ Eq and 



(0) 



Va,t3 



,(0) 
a,f3 



. Within RMT the quantities e^ 



and v^ L are 



uncorrelatcd. One might then conclude from Eq. H57(l 
that, in the limit u —^ 0, the gap distribution is a convo- 
lution of a Wigner distribution (evaluated with the GOE 
value of A without interaction) and a Gaussian whose 
width is w cr = \/2n 



1 cr(«i+i,n)V6w 



[there are 2n — 1 
interaction matrix elements in Eq. 157(1 ] . 

When compared with our previous results, the width 
a of the Gaussian would then be enhanced by the large 
factor \/2n — 1. This increase is not compensated by the 



FIG. 15: The peak spacing distribution -P(A2) (histogram) for 
the "local" interaction ensemble with u = 0.05 is compared 
with a convolution of a Wigner distribution with a Gaussian 
whose width is a{vl^^i „) (solid line). The dashed line de- 
scribes a convolution of a Wigner distribution and a Gaussian 



with a width of v2r; 



1 



(0) 



AO) 



"n + l,n)- 



Inset: the distribution 



hn,„ (histogram) is approximately fitted by a 



of ft, 

Gaussian (solid line) . The dashed line is the Wigner distribu- 
tion of En+i — e„. The number of electrons is n = 10. 



small decrease in the width of the Wigner distribution 
(the GOE value of A is a little smaller than the HF value 
Ahf)- The result is a distribution which is much wider 
than the one found in the numerical simulation. We show 
in Fig. El the gap distribution P(A2) (histogram) for a 
small value of u (u = 0.05). The solid line, predicted 
by Eq. (|75|l or Eq. H77I) . is a convolution of a Wigner 
distribution (with the appropriate A at u = 0.05) and a 

Gaussian with a width of cr = cr(w„_j_]^ „). As expected, 
the agreement is very good. However, the distribution 
that one infers using Eq. H57|) in the limit u ^ (dashed 
line) is much broader. The discrepancy is due to the 
fact that even for very small values of u, the difference 



,(0) 
''n+l.n+l' 



hnjn is strongly correlated with the HF matrix 
elements. Hence, the assumption made in using Eq. I|57|) 
is not correct. The inset of Fig.lT^lshows the distribution 



oih. 



(0) 
n+l,n+l 



hn,n (histograms). This distribution can 



be approximately fitted by a Gaussian (solid line) , but is 
distinctly different from the Wigner distribution of e„_|_i — 
e„ (dashed line). 



D. Peak— height distribution 

To measure the conductance, it is necessary to couple 
the dot to external leads. For weak coupling (i.e., for an 
almost isolated dot) and at low temperature (T ^ A), 
the conductance peak height can be expressed in terms of 
the ground-state wave functions of the dot with n and n+ 
1 electrons, and we can use the formalism developed here 
to study the statistical properties of the conductance. In 
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terms of the rates F' and V^ for an electron to tunnel 
into the dot from the left and right leads, respectively, 
the peak height is given by 



Goc 



T^/pr 



hkT r' + r*- 



n) 



The rates are squares of partial-width amplitudes, F'^'' = 
|^/,r|2^ In the CI model, the 7's are proportional to the 
projection of a single-particle eigenfunction (occupied by 
the electron that tunnels into the dot) on the respective 
lead (or point contact r). For symmetric leads that are 
separated by a distance large compared with the Fermi 
wavelength, the amplitudes 7' and 7*" are uncorrelated 
Gaussian random variables with zero mean value and a 
common second moment. This leads to the well-known 
predictions for the peak-height distributions^ 

A study of the peak height statistics in a small dot 
described by an Anderson model plus Coulomb interac- 
tions suggested that the peak-height distribution is only 
weakly affected by residual interactionsi^ To find out 
whether this result is generic, we use in the following the 
HF approximation to study the peak-height distribution 
within the framework of the induced two-body ensem- 
bles. 

In the presence of interactions (beyond charging en- 
ergy), we have 



7 ex ($(n-H)|V'Hi-)|*(ri)) , 



(82) 



where ip''{r) creates an electron at the point contact r, 
and $(n) is the ground-state wave function of the dot 
with n electrons. In the HF approximation, $(n) is a 
Slater determinant of the lowest n single-particle HF 
wave functions f/;! i ■ • ■ i '0n • 

In Koopmans' limit, the single-particle HF wave func- 
tions do not change with n. Expanding '0^(r) — 

"^x^K^^^x where a\ creates an electron in the single- 
particle HF state tpx of n electrons, we find 

7«C+i(r) (83) 

for the 71 ^ n + 1 transition. This expression is similar 
to the expression found for the CI model except that 
the single-particle HF wave function replaces the non- 
interacting wave function. 

We have shown in Section IIII El that the statistics of 
the single-particle HF wave functions are the same as 
in standard RMT. In particular, the distribution of the 
components ipn+iU) of the HF wave function ipn+i in a 
fixed basis \j) is given by Eq. (|65|l . For a chaotic bal- 
listic dot such a fixed basis Xi(i") is given by the free 
particle states with energy e = h'^k'^/2m, e.g., circu- 
lar waves Xj('') o^ Jj{kr)e^^^ {j = 0, ±1,±2...) where 
Jj are Bessel functions of the first kind. A single- 
particle eigenfunction of the non-interacting dot can be 
expanded V'a(i") = X]i V'a(i)Xi(i') and, for a chaotic 
dot, the components ipa{j) follow RMT statistics. Ex- 
panding the HF wave function ipn+i in the same basis 



tpn+i{r) = J2j V';(j)Xi(r), we can view Vn-n(r) as a pro- 
jection of ipn+i on a fixed vector whose components are 
Xj (r) • Because of the orthogonal (unitary) invariance of 
the HF ensemble, such a projection has the same distri- 
bution as the distribution (|66() of an eigenvector compo- 
nent. Thus the distribution of F — I7P is given by Ht)7|) 
and approaches a Porter-Thomas distribution for large 
m. Using ■0*_|_]^(i)-0„_(_i(j) = 5ii/m, we also find (see, for 
example. Section V.B in Ref. 



^i ryl 



a{Y)cr{Y) 



Jo(fc|i"i -i"r|) , 



(84) 



where a{'j) is the standard deviation of the partial am- 
plitude 7. Thus the correlation between the left and 
right partial width amplitudes decays in magnitude as 
~ (fc|r/ — r^l)"^/^, and can be ignored when the distance 
between the left and right leads is large compared with 
the Fermi wavelength. 

Since F' and r*" are uncorrelated and each follows a 
Porter-Thomas distribution, we conclude (in Koopmans' 
limit) that the peak-height distribution for the orthogo- 
nal (unitary) many-body ensembles of Section ^ is the 
same as the distribution derived from standard RMT. 

Beyond Koopmans' limit, the single-particle HF wave 
functions depend on the number of electrons on the dot 
and are accordingly denoted by i/Jx ■ Expanding in the 
HF basis of n electrons, ^p^r) — X^a^a ('')°'A (where 
a\ creates an electron in tpx ), we find 



/ . tX 

A=n+1 



($HF(n+l)|7^t(r)|ci,HF(n)) ^i-rY] vA"^*(r)detA(^), 

(85) 
where A^^) is an (n + 1) x (n + 1) matrix defined by 

4'j - (V^i"+'VJ"') and Ai^Ui = (^i"+^)|^i"') for 
a = l,...,n-t-l; $ — l,...,n. 

An equivalent expression for 7 can be derived by 
rewriting 7 ex (<I>hf(»^)|^(i")|^hf(«+ 1))* and expanding 
in the HF basis of 71 + 1 electrons tpir) — ^-^ ipx ('')^a 
(where now ax annihilates an electron in tp^^ ) . We find 



($HF(n+l)|Vnr)|*HF(n)) 



n+1 
t/„M* /..N\ V^, 



^(-)'-VA"+'^*(r)detBW, 



A=l 



^(86) 
where B^^^ is an n x n matrix defined by B — 



(V-i"^^^ I^AJ)"^) for a = 1, . . . , a - 1, a + 1, . . . , n + 1 and 
/3= l,...,n. 

Taking for 7/13^"' (r) a projection of the HF wave function 

-0)^" on a fixed vector, we have used Eq. H85|) to calcu- 
late the distribution of the renormalized partial width 
F = F/F in the HF approach (alternatively, we can use 
Eq. (|86|l with 7/13^" (r) given by a projection of the 
HF wave function i\r^ ). Fig. 1161 shows (symbols) the 
distribution P(lnf) versus Inf for u = 0.05,0.1,0.15 
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and 0.25. The solid line is the distribution (|67|l for 
TO = 30 and the dashed line is the Porter-Thomas distri- 
bution. For li = 0, only the first (last) term in Eq. H85(l 
(Eq. (ISHl) differs from zero, and the distribution is given 
by Eq. H67|l . As u increases, more terms contribute and 
the distribution of V gets closer to the Porter-Thomas 
distribution (this is plausible because of a central-limit 
theorem). The inset of Fig. Iltil shows the average value 
of r in units of the average value of P in the CI model. 
This average value decreases from 1 as u increases. We 
have also calculated correlations of partial widths (e.g., 

(fip^- f'f'^)/[cr(r'-)cr(P')]), and found them to be gen- 
erally weaker than predicted by Eq. (|69|l . 



0.3 







InP 



FIG. 16: The distribution P(lnf) versus Inf where f is the 
renormalized partial width in the HF approximation. Shown 
are results for u = 0.05 (triangles), 0.1 (squares), 0.15 (circles) 
and 0.25 (diamonds). We use the ensemble 1)17^ with m = 
30 and n = 10. The solid line is the distribution 16711 and 
the dashed line is the Porter-Thomas distribution. We have 
magnified a section of the graph to show more details. Inset: 
T/T{u — 0) versus u. 

We conclude that in the HF approach, the partial- 
width distribution is not affected by the fluctuations of 
the residual interaction and (in the limit of large m) 
remains Porter-Thomas like. Consequently, the peak- 
height distribution is similar to the distribution predicted 
from standard RMT. The residual interaction only affects 
the average values of the partial widths. 



V. SUMMARY AND CONCLUSION 

In this paper we have discussed a generic approach 
towards the understanding of the statistical properties of 
an (almost) closed diffusive or chaotic quantum dot in 
the limit of large but finite Thouless conductance g. 

The randomness of the single-particle Hamiltonian of 
a closed diffusive or chaotic dot induces randomness into 
the two-body interaction matrix elements when these are 



expressed in the eigenbasis of the single-particle Hamilto- 
nian. In the first part of this paper we have classified the 
resulting induced two-body ensembles. These depend on 
both, the underlying space-time symmetries of the dot, 
and the symmetries of the interaction matrix elements 
under permutations of the single-particle basis. Our clas- 
sification applies to the non-antisymmetrized matrix el- 
ements and therefore holds in the presence of spin de- 
grees of freedom. The ensembles for the spinless case 
follow directly by antisymmetrizing the interaction ma- 
trix elements. We have ignored spin-orbit scattering in 
the dot, so the resulting two-body ensembles have either 
orthogonal or unitary symmetry depending on whether 
time-reversal symmetry is conserved or broken. Aside 
from the symmetries, the ensembles are characterized by 
three parameters: The mean level spacing A due to the 
random one-body Hamiltonian, the parameter u which 
measures the fluctuations of the two-body interaction 
matrix elements, and the number n of electrons on the 
dot. 

The presence of two-body interactions poses difficul- 
ties in treating the many-body system. Therefore, in the 
second part of this work we have used the HF approxi- 
mation to study generic interaction effects in an (almost) 
closed dot. The HF approximation is not only tractable, 
but also allows us to use a single-particle formulation, 
optimized to include interaction effects. 

Applying the HF approximation to the two-body en- 
sembles, we have developed a generic statistical HF ap- 
proach. In this approach we have solved the HF equa- 
tions for a large number of realizations of the ensem- 
ble, thereby generating a one-body ensemble of single- 
particle HF energies and wave functions. We have stud- 
ied the statistical properties of this induced HF one-body 
ensemble, distinguishing the lowest n filled levels from the 
remaining empty ones. For both the filled and the empty 
levels separately, the nearest-neighbor level spacing fol- 
lows the Wigner distribution. This is not true for the HF 
gap separating the filled and empty levels. The gap dis- 
tribution has the form of a convolution of a Wigner dis- 
tribution with a Gaussian whose width is proportional to 
u, while the average value of the gap increases quadrat- 
ically with u (in leading order). By studying a suitable 
spacing correlator, we have found that the gap weakens 
the correlation between filled and empty levels, but simi- 
lar correlators within the filled or empty levels separately 
continue to follow RMT. We have shown that the single- 
particle HF ensemble satisfies orthogonal (unitary) in- 
variance in the presence (absence) of time-reversal sym- 
metry. Consequently the HF wave function components 
(in a fixed basis) satisfy RMT statistics. We have also 
studied the statistics of the interaction matrix elements 
in the HF basis (which reficct the statistical properties 
of the HF wave functions). The distributions of the HF 
matrix elements are Gaussian and, thus, can be char- 
acterized by their first two moments. The average of a 
diagonal interaction matrix element acquires a u^ correc- 
tion. The value of that correction depends on whether 
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the corresponding single-particle orbitals are both filled, 
both empty, or whether one is filled, the other, empty. 
The variances, on the other hand, remain close to their 
values in the non-interacting basis. 

The HF ensemble can be used to study the statistical 
properties of various observables. In particular, we have 
studied the generic properties of the peak-spacing dis- 
tribution and of the peak-height distribution. A simple 
interpretation of the results is provided in Koopmans' 
limit where it is assumed that the HF wave functions do 
not change upon the addition of an electron to the dot. 
The peak-spacing distribution is well approximated by 
a convolution of a Wigner distribution with a Gaussian 
whose width is the standard deviation of a diagonal inter- 
action matrix element. The peak-height distribution, on 
the other hand, is found to be insensitive to the residual 
interactions. 

We have confined ourselves to the HF approximation 
for spinless electrons. The presence of spin leads to ad- 
ditional technical difficulties in a mean-field approach. 
In particular, only the z-component of total spin (but 
not total spin itself) is conserved. Because of this re- 
striction in our calculations, it is not possible to make a 
quantitative comparison between our results and exper- 
iment. In the presence of spin, the single-particle lev- 
els are doubly degenerate. Therefore, the peak-spacing 
distribution is expected to be bimodal. The exchange 
interaction associated with the spin degrees of freedom 
explains the finite-temperature suppression of the width 
of the peak-spacing distributioufi^ but at low tempera- 
tures the distribution is expected to remain bimodal. No 
such bimodality is observed in the experiments, however. 
This fact is related to the fluctuating part of the inter- 



action (which is of order u oc A/g). We have seen that 
in the spinless case, these fluctuations affect the shape of 
the peak-spacing distribution and lead to a distribution 
which is intermediate between a Wigner distribution and 
a Gaussian distribution, in qualitative agreement with 
experiments. It would be interesting to study the corre- 
sponding finite-g effects in the presence of spin, using the 
HF approximation within the induced two-body ensem- 
bles discussed in Section|n] This should lead to a generic 
peak-spacing distribution that can be directly compared 
with the experimental data. As for the peak-height dis- 
tribution, RMT plus exchange interaction (without the 
fluctuating part of the interaction) were found to provide 
a quantitative description of the observed features at low 
temperaturesii^i^ We have shown here that in the ab- 
sence of spin, the peak-height distribution is insensitive 
to the fluctuating residual interaction. If this conclusion 
continues to hold in the presence of spin, it will fully ex- 
plain the observed low-temperature peak-height distri- 
bution. A definitive answer could be provided using the 
induced two-body ensembles of Section ^ that include 
spin. 
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